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Abstract 

Incorporating a realistic model for accretion of ultra-relativistic 
particles by primordial blackholes (PBHs), we study the evolution 
of an Einstein-de Sitter universe consisting of PBHs embedded in a 
thermal bath from the epoch ~ 10~ 33 sec to ~ 5 x 10 -9 sec. In this 
paper we use Barrow et al's ansatz to model blackhole evaporation in 
which the modified Hawking temperature goes to zero in the limit of 
the blackhole attaining a relic state with mass ~ m v \ . Both single mass 
PBH case as well as the case in which blackhole masses are distributed 
in the range 8 x 10 2 - 3 x 10 5 gm have been considered in our analysis. 
Blackholes with mass larger than ~ 10 5 gm appear to survive beyond 
the electroweak phase transition and, therefore, successfully manage 
to create baryon excess via X — X emissions, averting the baryon 
number wash-out due to sphalerons. In this scenario, we find that the 
contribution to the baryon-to-entropy ratio by PBHs of initial mass 
m is given by ~ e£(m/l gm) -1 , where e and £ are the CP-violating 
parameter and the initial mass fraction of the PBHs, respectively. For 
e larger than ~ 10~ 4 , the observed matter-antimatter asymmetry in 
the universe can be attributed to the evaporation of PBHs. 
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1 Introduction 



That the Milky Way is essentially made of matter is evident not only from 
the landings of space probes on Moon and other planets without any disas- 
trous consequences but also from the absence of anti-nuclei in the observed 
cosmic rays, and from the observations of Faraday rot at ion 0]. Observational 
support for absence of significant quantity of anti-matter beyond our Galaxy 
exists, but it is of indirect nature Since visible mass in the universe is 

chiefly in the form of baryonic matter, the inferred matter-antimatter asym- 
metry essentially boils down to the problem of the origin of baryon asym- 
metryQ The baryon asymmetry is characterized by the baryon-to-photon 
ratio r] = rig/ny, with and n 7 being the number densities of net baryons 
and photons, respectively. According to standard big-bang nucleosynthesis 
calculations, the predicted abundances of light elements depend only on the 
free parameter 77 and are in apparent agreement with the observed abun- 
dances provided rj lies in the range (2.8 — 4.5) x 1CT 10 ||. Recently, Tytler et 
alQ] have estimated the baryon-to-photon ratio from the observations of deu- 
terium abundance in a high red-shift quasar absorption system and according 
to their measurements, log 77 = —9.18 ± 0.4 ± 0.4 ± 0.2. 

The esthetically appealing scenario of the universe consisting of equal 
amount of baryons and anti-baryons at the instant of creation is still com- 
patible with a non-zero rj if one invokes Sakharov conditions, namely, of 
having B, C and CP violating interactions in out-of-thermodynamic equilib- 
rium condition sometime in the early history of the universe f|. The Grand 
Unified Theories (GUTs) of fundamental forces incorporate baryon number 
violating interactions naturally while CP violation can be introduced in such 
theories in many different ways (it is to be noted that CP- violation added 
theoretically in GUTs, in general, is not related to the observed CP-violation 
in the K° — K° system 0) and therefore it is not surprising that GUTs pro- 
vide a natural framework for the generation of baryon asymmetry through 
decay of X — X bosons J7[. However, through the work of Kuzmin, Rubakov 
and ShaposhnikovQ it came to be appreciated that baryon number violation 
can take place during electro-weak phase transition (EWPT) and that such 
processes could erase baryon-asymmetry generated prior to EWPT era. 

Use of 5-violation in electro-weak theories to produce excess baryons 

1 If neutrinos are massive then the gravitating mass may as well be dominated by 
leptons. However, there is hardly any direct measure of the lcpton number of the universe. 
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has also been made in the literature but one of the major obstacles in 
this scenario is the requirement of low Higgs mass which is in direct conflict 
with the experimental lower limit of m# > 88 GeV|TIJ. It appears that in 
the minimal version of electro-weak theory, generating baryon asymmetry 
may not be possible at all|l]J and especially with the discovery of the top- 
quark with a mass around 175 GeV[12| there is hardly any region left in 
the parameter space of the standard model to produce observed baryon-to- 



photon ratio p3 



The other major scenario of generating baryon asymmetry is to invoke 
Hawking evaporation of black-holes. The early sketchy ideas of Hawking 
and Zeldovich took proper shape with the advent of GUTs, giving rise to 
a picture of black-holes of small mass emitting X and X bosons thermally 
which subsequently decay and in the process violate B, C and CP, lead- 
ing to a production of baryon excess |n]]. At the fundamental level, this 
scenario has an attractive feature in that it combines ideas of black-hole 
thermodynamics [15] [16 1 on one hand and GUT on the other, to explain the 
observed matter-antimatter asymmetry in the universe. One of the impor- 
tant ingredients of this picture is the occurrence of mini-black-holes having 
mass less than 10 14 gm. It is obvious that such black-holes cannot emerge 



as end products of stellar evolution. However, Zeldovich and NovikovfT7 



and Hawking ||18|[ argued that primordial black- holes (PBHs) of small mass 
can be generated from the space-time curvature, and subsequently, CarrfT9| 



showed the possibility of creating PBHs from density fluctuations in the early 
universe. In the context of inflation, several authors have discussed mech- 
anisms to produce PBHs using the general idea that bubble wall collisions 
may trap pockets of false vacuum region that subsequently collapse to form 
black-holes ||20|| . In a recent work, Nagatani |21| has proposed an interesting 
blackhole-electroweak mechanism of baryogenesis that requires the presence 
of a blackhole to create a domain wall around it, leading to genesis of baryon 
excess without the need of a first order electroweak phase transition. 

Previous paragraphs of this section indicate that although GUTs can 
naturally generate baryon asymmetry, any baryon excess generated prior to 
electro-weak era is erased due to sphaleron transitions, while at the same 
time, creation of baryon asymmetry solely due to electro-weak processes is 
fraught with uncertainties as well as the requirement of low Higgs mass, con- 
trary to the experimental situation. Under the present circumstances, it is 
therefore natural to explore alternate means to explain matter-antimatter 
asymmetry. Since the existence of PBHs in the early universe is rather 
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generic, one ought to carefully re-examine the mechanism of generating baryon 
asymmetry through black-hole evaporation. In such a scenario, the crucial 
point to investigate is whether PBHs survive after the EWPT has taken 
place, so that the baryon asymmetry created due to their subsequent Hawk- 
ing evaporation survives, leaving an imprint till the present epoch. 

The present paper is an attempt to critically examine the evolution of the 
masses of a collection of PBHs created after the end of inflation, taking into 
account both the accretion of background matter by the black-holes as well 
as the mass loss due to Hawking emission. The paper has been organized 
in the following manner. In Section 2 we discuss processes responsible for 
the change in a black-hole's mass, and thereafter, we develop a formalism 
to describe accretion of relativistic matter by mini-black-holes. The subject 
of black-hole mass spectrum and its evolution is tackled next, in section 3, 
along with a discussion on the cosmological evolution of a mixture of PBHs 
and relativistic matter. Section 4 deals with the study of evolution equations 
numerically as well as a detailed analysis of the numerical solutions pertaining 
to the survival of PBHs past the EWPT. In section 5, we calculate baryon 
excess resulting from the decay of X — X bosons emitted by the PBHs during 
their final stages of Hawking evaporation, and then discuss the implications 
of these results to the question of matter-antimatter asymmetry. Finally, we 
end with a brief discussion of the above scenario in section 6. 



2 Evolution of the mass of a black-hole 

2.1 Mass loss due to evaporation 



Bekenstein's conjecture[|T3J that the area of the event horizon of a black-hole 
being a measure of its entropy was vindicated by the classic work of Hawk- 
ing in the early seventies who showed that when quantum effects around a 
black-hole are included, the black-hole emits particles with a thermal dis- 
tribution corresponding to a temperature Tbh that is proportional to the 
surface gravity at the event horizon| I3|||16|| , and is given by the relation, 



m pl c 



2 



T BH = ^ T (1) 

where m and m v \ are the mass of the black-hole and the Planck mass, re- 
spectively. According to eq. (d), PBHs created in the early universe with a 
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mass m 10 14 gm would be decaying today in a burst of high energy radia- 
tion, and there exists in the literature, upper bounds on the abundance of 
such PBHs from the observed level of cosmic 7-ray flux. As pointed out by 
Zeldovich and others, the expression in eq. ([]]) for the Hawking temperature 
can only be an approximation and is amenable to modifications at Planck 
scale because of the effects of quantum gravity. In fact, particle physicists 
have shown from various angles that Hawking evaporation may cease when 
the black-hole reaches the Planck mass scale leading to a massive relic. In 
this context, an interesting toy model inspired by superstring theories has 
been considered by Barrow et al|22[ in which the expression for the black-hole 



temperature has been modified by including correction terms that contain 
powers of black-hole mass in units of Planck mass. Following Barrow et al's 
ansatz, one can therefore express the black-hole temperature as 



rp _ "V 

J-BH — — — 



57T 



m \ m J 



,,2 



¥ < 2 > 



where k is a non-negative constant. For n > 2 and k ~ 0(1), it is clear that 
for holes of mass m ^> m p i, eq. ([!]) is a limiting case of eq. (|2|). According 
to eq. (0), as the hole mass decreases due to evaporation, initially there is 
a rise in the hole's temperature but as m approaches m p i the temperature 
starts falling and becomes zero when the mass of the hole reaches the value 
m re/ = ^^ n ~^vn v i. Therefore, Barrow et al's ansatz implies stable black-hole 
relics of mass m re i ~ m p i. To estimate the rate of mass-loss from eq. (0), we 
may work in the frame-work of radiative transfer, assuming that the hole's 
event horizon acts like a perfect black-body surface. In such it is 

easy to show that the energy flux F is related to the energy density e at the 



surface of interest in the following manner |23 

The effective energy- density of ultra-relativistic particles due to Hawking 
evaporation in the vicinity of the event horizon is related to the temperature 
of the black-hole by 

2 BH 7,4 

where gf = g^ H + (7/8)gf H is the effective number of degrees of freedom 
at the temperature T BH , and g^ H and gf H are the corresponding degrees of 
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freedom for bosons and fermions, respectively. Therefore, the rate of mass- 
loss from the event-horizon is given by 



dm 
~dt 



-«2m 



4nR 2 s 



rripi 
m 



— K 



rripi 
rn 



(5) 
(6) 



where «2 = g* H c 2 /(30720irh) and k ~ 1. In arriving at eq. (|6]), we have 
made use of eqs. as well as the standard result for the Schwarzschild 

radius Rs = 2Gm/c 2 . The calculations that led to eq.(||) were based on 
modeling the black-hole event-horizon to be the surface of a black-body of 
radius Rs at a thermodynamic temperature Tbh- It is> therefore, interesting 
to compare our result with that of Don Page[^] which is based on rigorous 
numerical computations for black-holes of mass m > 10 17 gm. According to 
his calculations, the mass-loss rate for such holes is 



dm 
~~dt 



-2.011 x 10" 



he 4 



G 2 m? 



(7) 



If we assume that eq.(|7]) is valid also for 10~ 2 gm < m < 10 17 gm, then 
comparing (§) and (|7|) one obtains g+ H ~ 20, which is not too unreasonable 
since for holes of mass 10~ 2 gm one expects g+ H to be as high as ~ 100 (in 
most GUTs). 



2.2 Accretion of relativistic matter by a mini black- 
hole 

The temperature of the universe is expected to be extremely high just after 
the end of inflation, and therefore matter during that period will be in the 
form of ultra-relativistic particles. For particles with de Broglie wavelength 
A Rs, the capture cross-section corresponding to a Schwarzschild black- 
hole is ~ nr 2 where r c = (3-y/3/2)-Rg [p5| . When the de Broglie wavelength of 
a particle is larger than R s , the capture cross-section is likely to be negligible 
as the blackhole sees an incident wave rather than a point particle. For high 
energy particles with A -C Rs, we will make use of the geometric optics ap- 
proximation in which any such ultra-relativistic particle hitting a fictitious 
sphere of radius r c around the hole will be absorbed. 

If I v represents the specific intensity of such particles corresponding to 
energy hv and if dA is an area element on this fictitious sphere then the rate 
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at which energy is accreted by the hole per unit range of v per unit area is 
given by 



fdttcos6I u =ttI u (8) 



dtdvdA 

Since the effective area of capture is 47rr^, the rate at which energy is accreted 
in the frequency range [u, v + dv\ is given by 

^ = (2nr c ) 2 I u du (9) 

To obtain the total rate of accretion of energy we integrate eq.(||) over 
frequency keeping in mind that geometric optics approximation requires the 
lower limit of integration v m i n to be a few times c/r c . For ultra-relativistic 
particles, momentum is p ~ hvjc so that the number density of particles of 
species A in the frequency-range {y, v + dv) takes the form 

4vr# A v 2 dv 

n A {y)dv = ^ ehv/kT±1 (10) 

where qa is the spin-degeneracy factor for the A th species and the +(— ) sign 
refers to fermions (bosons). In eq. ( |TUD T is the temperature of the universe. 
Therefore, the specific intensity I uA corresponding to the species A is given 
byfH 



j _ chvn A (v) _ 9A hu 

{ vA — 47r — ^ e hu/kT _|_ l 



(11) 



Making use of eqs.@ and (]TT|) , we can express the net rate of energy accretion 
by a hole in the following manner 



dE _ r2nr, 
~dt ~ 



[12) 



where z/ min = a.\cjr c is the lower frequency cut-off, ot\ being a number of 
the order of 10 (this takes care of the fact that only particles with A <C Rs 
are considered to have been captured by the blackhole). In eq. (|12|) , g^and 
gy m are the total bosonic and fermionic degrees of freedom, respectively, for 
the cosmic soup. These are to be distinguished from g^ H and gf H introduced 
in section (|2.1|) . 
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The rate at which the hole's mass grows as a result of accretion is 



dm 



dt 



7T 3 C 5 
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oo 



dx + 




e x + 1 



dx 



(13) 



where x min = hv min jkT. 

In obtaining the above equation, we have made use of a change of vari- 
able in eq.(|l2D along with r c = (3y/3/2) Rs. We note that Er appear- 
ing in eq.(fi"3|) is the energy- density of the background relativistic particles, 
sr = 7T 2 g™ ni (kT) 4: / (30h 3 c 3 ) , g™ m being the temperature-dependent effective 
spin-degeneracy factor and is equal to g^ m + 7/8gJ m . From eq.flTB]) it is 
evident that accretion plays an important role for massive PBHs at early 
epochs when the temperature of the universe is very high so that energy 
density Er of the relativistic particles is large while x min is small. This is 
easy to understand from a physical point of view in the sense that only when 
the temperature is large that there are sufficient number of particles with 
de Broglie wavelength much less than the Schwarzschild radius of the PBHs, 
ready to be accreted. By the same token, when the hole-mass reaches a size of 
the order of m p i, neither accretion nor quantum evaporation is significant. 

3 Blackhole mass spectrum and evolution of 
the universe 

It is evident that the mass distribution of PBHs is intimately linked to the 
mechanism of their production. Several authors [[LJ], |2(], |26|, in the lit- 
erature have discussed blackhole mass spectrum from diverse angles. Since 
the mass spectrum is sensitive to production mechanisms and, since so far 
no particular model of PBH creation has been singled out, we adopt a very 
general procedure in this paper to analyse the evolution of blackhole mass 
spectrum. 

We consider a distribution function N(m,t) such that N(m,t)dm repre- 
sents number-density of PBHs with mass in the range (m, m + dm) at the 
cosmic epoch t. We assume that the creation of PBHs stopped after a cosmic 
epoch tpbt so that at later times in a given comoving volume the number of 
holes remain the same while their masses change due to a combination of 
Hawking radiation and accretion of background matter. Note that we are 
working under the assumption that the ultimate state of a PBH along the 
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course of its evolution is a stable relic of mass ~ m p i, i.e. a hole does not 
disappear completely as the original Hawking radiation mechanism would 
demand. Also, since the mass m of a hole changes with time, the mass 
distribution function at time t and at time t + dt are related as 

a 3 (t)N(m, t)dm = a 3 (t + dt)N(m, t + dt)dm! (14) 

where ml is related to m through ml = m + mdt and a(t) is the FRW scale- 
factor at cosmic epoch t. Making a Taylor expansion of quantities at the 



RHS of eq.(|l4]), and using the relation 



COtti \ 
1 + — dt (15) 
dm J 

we obtain 

dN ad,, , _ 

With the help of the mass distribution function N(m,t), we can also obtain 
an expression for the mass-density associated with PBHs as 



pBiiit) = \ mN(m,t)dm (17) 

Jm re i 

It is useful to express the black-hole mass-distribution as 

N(m,t) = N (t)f(m,t) (18) 

where Nq oc a~ 3 (t) so that J™ f(m,t)dm is independent of time. With the 
help of eq.(^) it can be easily shown that eq. ([l6|) reduces to 

Essentially, f(m,t)dm represents the number of black-holes with mass in 
the interval (m, m + dm) in a unit coordinate volume at the cosmic epoch 
t, while the dilution of black- hole number density due to the expansion of 
the universe is taken care of by the factor N (t) = A/a 3 (t). Differentiating 
eq.(|T7D with respect to t and then making use of eqs.(^) and (|T9| ) it can be 
shown that 

dpBH , 3a 
~dT 



3d f°° 
+ —Pbh = N mf{m,t)dm (20) 

CL J m re i 
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Since the total energy-momentum tensor is divergence free, we also have the 



equation 28 



2 d 
° dt 



(Pr + Pbh)ci 3 + 3p R a 2 d = (21) 



where pR = Er/c 2 is the mass density of radiation. Here we have assumed 
that the black-holes possess negligible peculiar speeds so that their contri- 
bution to pressure is insignificant. Using pr = c 2 pr/3 and eq.(|2T|) in eq.(pO|) 
we obtain: 

dot) d r°° 
1 R + A-pr = -N / mf{m, t)dm (22) 



dt 



Eq(|22j) just reflects, as is to be expected, the fact that an effective black-hole 
mass loss (or gain) would imply pr oc a~ 4 ~ a where a(t) is negative (positive) 
because of black-holes acting as source (sink) of radiation. 

In any mechanism of PBH production, the actual masses of the black- 
holes will be distributed in a discrete fashion, and therefore without loss of 
generality the distribution function can be expressed as 



K 



f(m,t) = 1 £l3iS{m-m i (t)) (23) 



8=1 



where fa are constant weights corresponding to rrii, and K is the number of 
distinct black-hole masses. It can easily be ascertained that the distribution 
function in eq(|23"l) indeed is a solution of eq( [L9|) , since 

d TYl ' 

—5{m - rrnit)) = ~ ™i(t)) (24) 

at m — rrii(t) 

and 

3 ^ rfi ■ 

dm r-J m — rrii{t) 



Consequently, we have: 



' il>R + A P R- = -N a {t)j2fam l (26) 



dt a ; ! 



The manner in which the individual mass rrii of a PBH changes with time 
depends on the combination of Hawking evaporation rate and the accretion of 
background relativistic matter as discussed in section pi. Therefore, making 
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use of eq and eq 
following result: 



in the context of a PBH with mass mjwe get the 



drrii 



2^2 



Qb'"" f 00 x 3 dx _i_ 9f roo x 3 dx 

~uni I 7> ,-,x 1 ~~ T~ n !ini , / t - _i_1 



lx_ , roo 



(27) 



From eqs.(17),(18) and (23), the mass density pbh associated with the PBHs 
can be written as: 



A 



p BH {t) = iV„(t)^ft mi (t) 



(28) 



The evolution of the scale-factor a(t) then follows from the flat FRW Einstein 
equation: 



-)' 

aj 



8ttG 



A' 



PR + N (t)J2(3 imi (t) 



i=i 



(29) 



In writing down the above equation, we have adopted the inflationary paradigm 
according to which universe in the post-inflationary phase is described essen- 
tially by a flat FRW model. In this paper, the evolution of the universe is 
determined by three coupled differential equations fl26|) ,(27) and ( ^9|) along 
with the fact that N (t) oc a~ 3 (t). 



4 Numerical Evolution 



In this section, we solve 2 + K coupled non-linear, first order differential 
equations ( p6f) , (27) and (p9| ), set up in the preceding section, numerically 
using Hemming's fourth-order, double precision predictor-corrector method. 
To begin with, we fix N (t) by demanding that j3imi(to)No(to) represents the 
initial fraction Q of total mass density p(to) that lies in blackholes having 
initial mass mi(t ) so that, 



Pimi(to)N (t ) = CiP(*o) 



(30) 



As N (t) oc a 3 (t), we have, from eq((3 



N (t) 



a 3 (t ) CiP(to) 
a 3 (t) Pirrii(to) 



(31) 
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Since N Q (t) is independent of i, we obtain the following relation between Q 
and fa, 

Ci oc fami(t ) (32) 



where the constant of proportionality in eqfl3"2]) can be determined from the 
following identity, 



K 



pn(to) = p(t ) -J2(ip(t ) 
leading to the following expression, 



(33) 



i=l 



K 



proportionality constant = y^ffim^tc 



\i=l 



PR^o) 
P(t ) 



(34) 



Substituting eq(|31"D in eqs fl26|) and (^), we obtain, 



dpR , A a 
—7- + ^Pr- 
at a 



8nG 



PR 



a 3 (to] 

a 3 (to, 
' a 3 (t) 



K 



p(*o)Eo 



VTLi 



i=l 



mi(t c 



■ p(to) S°^)j 



(35) 
(36) 



respectively. 

In our original formulation (see section 3), the blackhole initial mass- 
spectrum was completely specified by the set of numbers {fa, rrii(t ); i = 
1, . . . K}. Equivalently, since fa and Q are related by equation (|32]) , we may 
as well specify the spectrum by the set {Q, mj(t ); i = 1, . . . K}. 

For the purpose of numerical evolution, it is convenient to cast equations 



(37) 
(38) 

(39) 

(40) 



(27) , (pq) and ( j36"D in terms of dimensionless quantities defined below 

T 

a(r) 
R(t) 



t\jGpo 
a(t) 
a 

P R (t) 



-a 



MAt) 



Po 
rrii(t) 

(t ) 
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where p = p(to) an d °o = a(*o)- I* 1 terms of the above quantities, the system 
of differential equations assumes the following form: 



a = — 



R' 
M' 



K 



\ 6 i=i 



K 



-a £ c,m; 



mi{t Q )M?{Gp Q ) 



-1/2 



405 

7T 3 C 3 



G 2 p$Ra Ji — «2 



mi(to) ^ 



— K 



rripi 



mi(to) Mi 



(41) 

(42) 
(43) 



where prime denotes differentiation with respect to r and for convenience we 
have introduced 



Ji = J (x Qi , T) 



rOC 

Jxo, 



00 x 3 dx 



9b 



n uni I p x 1 

9+ Jx o, e — i 
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x 0i 
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SVZGrriiit) 
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00 x 3 dx 



e x + 1 



We choose to to be the cosmic-epoch when inflation ends ~ 10~ 33 sec, and set 
p = 10 56 GeV 4 , which is the density expected at GUT scale. In our numer- 
ical evolution program, the actual values used for the following parameters 
are listed below: 





= 10 
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= 0.1 
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= 3 




- n BH 
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= 50 
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First we consider the case when K — 1, i.e. at the end of inflation a 
fraction ( of matter lies in blackholes, all with initial mass m = m (to)- We 
study different models by varying ( in the range 10~ 3 to 10 _1 while mo runs 
through the range 10 3 to 5 x 10 5 gm. In fig (1) we plot a(t) for a typical choice 
of C and m . The plots of a(t) for a radiation-dominated (RD) FRW universe 
(a ~ t 1 / 2 ) and a matter-dominated (MD) FRW universe (a ~ t 2 / 3 ) are also 
given in the same figure. The initial behaviour of the system is that of a RD 
universe but soon the evolution of the scale-factor a becomes similar to that 
in an MD universe, and subsequently, as the PBHs evaporate, the dynamics 
becomes RD again. This is because initially energy density of relativistic 
matter gets depleted owing to accretion by PBHs, resulting in its decrease 
faster than the kinematic rate a -4 (see fig (1)) so that the dominant contri- 
bution from 'dust' like blackholes drives a faster expansion rate. We have 
also compared our results with approximate estimates obtained by assuming 
that a(t) ~ t 2/3 from t = t (end of inflation) to t = t E wpr (epoch of EWPT) 
and that a(t) ~ t 1 / 2 afterwards. For the range of parameters 10 3 < m < 10 5 
(gm) and 0.001 < ( < 0.1, the estimates agree with our numerical results to 
within an order of magnitude. 

In fig (2) we plot a typical mass m as a function of time. It is evident 
from the figure that growth of blackhole mass due to accretion takes place 
only in the initial period when the temperature and density of the universe 
is very high. This is anyway expected since the de Broglie wavelength A of a 
typical particle just after the end of inflation is ~ 10~ 28 cm, while the Rs for 
a blackhole of mass as low as ~ 100 gm is ~ 10~ 26 cm leading to a substantial 
accretion because of A < Rs criteria. At intermediate times the curve flattens 
out reflecting a balance between accretion and Hawking evaporation. During 
this phase, the dynamics is essentially MD since radiation loses out in the 
competition because of the expansion of the universe as well as its attenuation 
due to accretion by PBHs. 

Towards the end, Hawking evaporation begins to dominate the evolution 
of PBH mass as the accretion automatically gets switched off due to the de- 
crease in temperature and density of background radiation. In fig (3) we plot, 
for a typical choice of parameters ( = 0.01 and m = 2.5 x 10 5 , the ambient 
temperature of the universe T as well as the Hawking temperature Tbh of the 
black- hole. The straight line portion of the curve has slope equal to —2/3. 
Thus T falls, at intermediate times, as if the dynamics of the universe was 
akin to that of a MD universe. At later times, when evaporation becomes 
the dominant process in the evolution of the holes, the universe at first starts 



14 



cooling at a slower rate, but eventually reheats due to rapid evaporation of 
the blackholes (the reheat portion is not included in the figure). The point 
at which EWPT occurs is marked by an arrow in the figure, (i.e. T is ~ 100 
GeV at this instant of time) corresponding to a value of ~ 2 x 10~ 13 sec. We 
note that the epoch of EWPT is considerably lower than the standard value 
of ~ 10~ 10 sec obtained from the time-temperature relation in big-bang mod- 
els. The reason for this is not hard to understand, as depletion of radiation 
by the accreting PBHs leads to a MD phase causing the T to decline faster 
than the usual t -1 / 2 fall. 

Now, the amount of reheating should be such that the temperature of the 
universe does not rise above the EWPT temperature (~ 100 GeV), because, 
otherwise the sphaleron processes will be re-ignited, leading once again to a 
washing out of BAU generated. This, in effect, constrains our parameters 
( and itlq. In fig (4) we plot the combination of ( and mo for which the 
re-heat temperature is 100 GeV, and these points are empirically fitted with 
a curve. From the numerical evolution, we find that the region lying below 
to the right of the curve consists of those values of (C^o) f° r which re- 
heat temperature remains below 100 GeV. While the region lying left of 
the curve consists of those combinations for which PBHs evaporate away, 
reaching the relic state before EWPT, and therefore are not of any use as far 
as baryogenesis is concerned. From fig (4) it is evident that for ( lying in the 
interval (10" 5 ,0.1) a PBH with initial mass less than ~ 2x10 gm converges to 
the relic state before EWPT, and hence does not contribute to generation of 
baryons. Therefore, we find that if the initial PBH mass spectrum is a delta- 
function peaking at the mass mo, baryogenesis through blackhole evaporation 
is viable only when the initial mass of the PBHs exceeds ~ 2 x 10 5 gm for 
reasonably low values of (. 

Next, we consider the case in which blackhole masses at time to are dis- 
tributed in a pseudo-Maxwellian manner as shown in fig (5). The blackhole 
masses fall in a range from 8 x 10 2 gm to 3 x 10 5 gm, with 22 distinct mass 
values contributing to a total fraction Y% 2 ~ 0.09 of the mass density of the 
universe just after the end of inflation. From the plot fig (6) it is apparent 
that blackholes with larger initial mass accrete background hot matter at 
higher rates than those with smaller initial mass, as expected from the fact 
that higher mass PBHs have larger cross-section for absorbing matter. We 
find that those PBHs with initial mass greater than 1.5 x 10 5 gm reach X — X 
emitting phase after the epoch 1.9 x 10~ n sec, the instant at which EWPT 
takes place for this spectrum of masses. Once again we find that EWPT oc- 
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curs sooner than that in the standard model. There are 7 of such blackhole 
masses which finally contribute to the production of baryon excess. Because 
of the wide distribution of blackhole masses, the instants at which the PBHs 
reach the relic mass are staggered, hence no sharp re-heating takes place in 
our analysis, rather the temperature of the universe falls at a slower rate till 
the largest size blackhole (with initial mass = 3 x 10 5 gm) evaporates, leaving 
behind a relic mass around the epoch ~ 3 x sec, when the temperature 
of the universe is ~ 9 GeV. The decline of temperature with time is shown 
in fig (7). 

Even in the case of blackhole mass distribution, K being larger than 
1, in principle, one can constrain the parameter space (G> ^(^o)) from the 
requirement of re-heating less than 100 GeV (as undertaken when K = 1, see 
fig (4)), however the exercise is enormously time consuming, and is beyond 
the scope of the present paper. 



5 Baryogenesis 



We saw in the previous section that for ( « 0.01, blackholes created with 
mass less than k 2 x 10 5 gm evaporate and reach the relic state before the 
EWPT and hence their contribution to baryon asymmetry is doubtful due to 
the expected 5-violation induced by sphalerons. However, blackholes with 
initial mass larger than « 2.5 x 10 5 gm certainly ought to be considered as 
sources of baryogenesis since they reach the X — X emission phase well past 
the EWPT. In this section, we proceed to estimate the quantity of excess 
baryons resulting from blackholes whose Hawking temperature reaches GUT 
scale after EWPT. 

Representing the specific intensity of X-bosons radiated with energy hu 
from a blackhole by I„ , we have the relation (e.g., see p3| ) 



rX 



<(n) 



(44) 



where u*(Q) is the specific energy density and v is the speed of the emanating 
X-bosons. 
With 

V = C 1 



mc 



hv 



(45) 



16 



and 



v c 4 e hu/kT BH (m) _ 1 



we may express I* as 

his 3 



I 



x 



v 2 
C z 



2\ 2 



hp 



^hi//kT B H(m) 



<A - (47) 



It is to be noted that gx and T BH {m) are the spin degeneracy factor of 
X-bosons and Hawking temperature of a blackhole of mass m, respectively. 
The flux-density of X-bosons at a distance r from the blackhole is given 

by 

Ff = (48 ) 

Therefore, from equations ( fi?|) and (|48| ) the rate of emission of X-bosons 
from a blackhole of mass m is derived to be 



dN x (m) r°° F* ■ ^r 2 dv 



dt Jm x c 2 /h hv 

4irR 2 s c 4 g x m x 



(49) 



h 3 

where m>x is the mass of the X-boson and 



HVi) (so) 



r°° v 2 — 1 

'<*> = / ^-jdv (si) 

while 

Since, at any given cosmic epoch t, the number density of blackholes with 
mass lying in the interval (m, m + dm) is N (t)f(m, t) (see eq (|18|)), the rate 
at which X and X bosons are generated in a unit proper volume is given by 

dn X x n „ u f dNxijn) 



dt 



2N (t)J^^f(m,t)dm (53) 



In eq fl53"D the factor 2 arises because we have included production of X- 
bosons as well. Making use of the form given in eq (|23|) we can express eq 
(|53|) as 

^f = 2«»^ (54) 
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The lifetime of a X-boson tx = T^ 1 turns out to be ~ 10 36 sec when nix ~ 
10 14 GeV |29| which is negligible in comparison with the time scales over 



which blackhole mass changes or the universe expands appreciably. Hence, 
the rate of increase of net baryon number in a unit proper volume is 



dn x x 
~ e 



dt 
with 

e s r(x -> g p - r(x -> g) 

being the net baryon number generated by the decay of a pair of X and X 



If riB(t) represents net baryon number density at the cosmic epoch t then 
j t {a\t)n B {t)) = ea\t) d -^ (56) 



Employing eq ( |3lD and (|53|) in ( pip and then integrating the latter, we obtain 
a 3 {t)n B (t) - a 3 (t EW p T )n B (t EWPT ) = 2e{p al) jl \ f ~~377~ L dt> 

i=l m i{t0) Jt-EWPT at 

(57) 

Assuming that prior to blackhole baryogenesis, the net baryon number in 
the universe is zero (i.e. nsitEWPT) = 0) and making use of eq (|50"D in eq 
([571) ) we get the following expression for the net baryon number density at 
any time, 

After EWPT has taken place, the evolution of PBH mass is totally dominated 
by eq (^) since the de Broglie wavelength A of a typical particle is larger than 
~ 10 -15 cm, while the Rs corresponding to a blackhole of mass as high as 
~ 10 7 gm is only ~ 10~ 21 cm. Hence, using eq (^) we can change the variable 
of integration in eq(j58|) from t' to m^t') so that, 

/ m 2 t (t')I( yi (t'))dt' = pl HimSowe^^iitEWPT)) (59) 

JtEWPT a 2 (87rmx) 4 
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where H is defined to be, 



H(m,i(lower), m» {t E wpr) ) 



i(tEWPT) 1 
^2 



e /yi ^ e 



fc=l K k=l K 



dm 



Imi(lower) ?/j 

(60) 

In obtaining eqs (|59"D and (|60[) we have used the series equivalent of the 
integral given in eq (|5l"l) .The value of m^lower) is set by requiring £/, to be 
10 -3 since the series given in eq(|60|) is negligibly small for smaller values of y^. 
This automatically takes into account the fact that only those PBHs matter 
for BAU that are capable of emitting X — X after EWPT. For PBH masses 
larger than 10 5 m p z, the value of H is 2.7 x 10~ 2 and becomes insensitive to 
the exact value of vn^EWPT) thereafter. Therefore, for the 7 PBHs that 
survive the EWPT, we have, 



22 j- 

E -^jrrHimiilower), rm(t EWPT )) = 1.97 x 1(T 9 (61) 

i=16 "WO J 

The entropy density of the universe at any epoch t is given by, 

• = £*- , 0& I *M (62) 

We estimate the baryon-to-entropy ratio at t ~ 3 x 10~ 9 sec, when all the 
PBHs settle on to the relic state, by making use of eqs (|58|) , (|59"D , (j60|) , fl6"l"|) and 

Mt) - ™ x io-»e 9x rev rev <«» 



s(t) * \ 100 J \100 

The contribution to baryon-to-entropy ratio by PBHs with initial mass mo 
and initial mass fraction ( goes roughly as, 

^^Jj-VVsq-Vsry 1 (64) 

s sy \igm) ^ ioo y v ioo y v ; 

Hence, in the case of a delta-function mass spectrum with £ « 0.01 and 
mo ~ 2.5 x 10 5 gm, one obtains a baryon-to-entropy ratio of ~ 4 x 10~ 8 e, 
with gx — 1. Thus one may use eq(|64"D along with the value oins/s ~ 10~ n , 
that follows from observations, to put a constraint on eQ/m^. This implies 
that one requires the CP-violating parameter e to be around ~ 10 -4 to 
generate excess baryons from evaporating PBHs. 
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6 Discussions 



To study the evolution of PBHs, in the early universe, that undergo accretion 
along with steady mass loss due to Hawking evaporation, we have laid down 
a formalism which can handle any blackhole mass spectrum that can be de- 
composed as a sum of weighted 5-functions. Accretion of ambient hot matter 
by a blackhole has been modeled in the limit of geometric approximation, so 
that only those particles with de Broglie wavelength less than about a tenth 
of Schwarzschild radius are considered for absorption by the blackhole. The 
evolution of a flat FRW universe and the PBHs has been studied numerically 
to find conditions under which blackholes survive past the electroweak phase 
transition in order that their subsequent evaporation leads to baryogenesis. 

The basic picture which emerges is the following. In the case of a blackhole 
mass spectrum that peaks sharply at a single mass value mo, when ( (the 
initial mass fraction of PBHs) is of the order of 1%, PBHs with initial mass mo 
less than about 2.3 x 10 5 gm evaporate before EWPT. Therefore, only PBHs 
with m greater than this critical value need be considered for generation 
of BAU. Here, we wish to point out that the model of accretion which one 
considers can make an immense difference in the final result of the analysis. 
If one uses a simple spherical model of accretion in which the capture-cross 
section is just 7ri?| and with no de Broglie wavelength based cutoff then 
blackholes of initial mass mo ~ 10 3 gm can successfully live past the EWPT, 
and eventually contribute to the BAU (see Majumdar et al. in [|14]]). While 
on using the same set of parameters with a wavelength based cutoff model 
of accretion, we find that PBHs of such small initial mass do not survive 
beyond the EWPT. 

For reasonable choice of parameters, we find that in the case of PBHs 
with a distribution of mass ranging from 8 x 10 2 - 3 x 10 5 gm, blackholes 
with initial mass larger than about ~ 10 5 gm reach the relic state much after 
EWPT. Because of the presence of blackholes with mass less than 10 5 gm that 
evaporate at a faster rate, pumping in energetic particles into the surrounding 
medium, the ambient temperature in this case declines at a slower rate, and 
hence EWPT takes place later than in the case when all PBHs had the same 
mass of 2.5 x 10 5 gm. As described in sections (^|) and ([|), the evolution 
of mass spectrum is totally determined by the manner in which individual 
blackhole masses change with time, /3j or equivalently Q remaining fixed for 
all times. As an illustration, we have shown the evolution of mass spectrum 
in fig (6) for a particular set of Q. 
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We wish to point out that accretion is important only during the initial 
stages just after the end of inflation when the temperature of the universe 
is ~ 10 13 GeV, causing an increase in the mass of a blackhole by a factor 
of ~ 4. There are two factors responsible for a blackhole of initial mass of 
~ 10 5 gm to live after the EWPT. One being the increase in the mass due 
to accretion, while the other is the occurance of EWPT sooner than that in 
a model in which there is no depletion of radiation due to PBHs acting as 
sinks. For blackholes with mass less than ~ 10 5 gm, accretion is less due to 
the reduction in capture cross-section because of which the rate of depletion 
of radiation is not large leading to a delayed occurance of EWPT, after the 
blackholes have reached the final relic state. 

Barrow et al's ansatzfE^] which has been used in this paper to take into 
account expected modification of Hawking emission becomes important only 
when the mass of evaporating blackholes fall below ~ 10m p i. Our numerical 
results are not sensitive to the exact form of modified blackhole temperature. 
For baryogenesis, significant quantity of X — X are emitted only during the 
phase when blackhole temperature is ~ Tout, because of which the integral 
H is not sensitive to the upper limit m(t EWPT ) so long as the latter is larger 
than 10 5 rripi. Therefore, the final expression for baryon-to-entropy ratio 
turns out to be rather simple (see eq (@)), implying that if at the end of 
inflation 1% of total matter goes into creating PBHs with initial mass 2.5 x 10 5 
gm then this scenario can successfully lead to BAU provided the CP-violating 
parameter e is over 10~ 4 . Thus production of baryon excess through blackhole 
evaporation is a viable alternative to GUTs or electroweak baryogenesis, 
although there is no denying that because of the presence of parameters like 
Ci and (to) whose values a priori are uncertain, this scenario cannot provide 
meaningful constraint on the value of e. 
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Figure Captions 



1. Figure 1: The evolution of the scale-factor a(t) for ( = 0.01, m = 2.5 x 
10 5 gm. The plots t 1 ' 2 and are provided for comparision. 

2. Figure 2: The evolution of the mass m(t) of the PBHs for a typical 
choice ( = 0.01, m = 2.5 x 10 5 gm. 

3. Figure 3: The temperature T of the background thermal bath and the 
Hawking temperature T BH for a typical choice ( = 0.01 and m Q = 
2.5 x 10 5 gm. The instant of EWPT is marked by an arrow. 

4. Figure 4: The combinations ( and mo for which the reheat temperature 
= Tewpt =100 GeV. The region with acceptable reheat temperatures 
< 100 GeV is indicated in the figure. The analytical fit with dotted 
line is purely empirical. 

5. Figure 5: Black-hole mass spectrum: plot of Q against mj(t ). 

6. Figure 6: The evolution of the masses m,i(t) of a collection of PBH 
masses distributed according to the spectrum shown in fig (5). 

7. Figure 7: The cooling of the universe for the case where PBH masses are 
distributed according to the spectrum displayed in fig (5). The epoch 
of EWPT is marked by an arrow, and it takes place at 1.9 x 10~ n sec. 
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Figure 1: The evolution of the scale-factor a(t) for ( = 0.01, m = 2.5 x 10 5 
gm. The plots a ~ t 1 / 2 and a ~ t 2 / 3 are provided for comparision. 
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Figure 3: The temperature T of the background thermal bath and the Hawk- 
ing temperature Tbh for a typical choice ( = 0.01 and mo = 2.5 x 10 5 gm. 
The instant of EWPT is marked by an arrow. 
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Figure 6: The evolution of the masses rrii(t) of a collection of PBH masses 
distributed according to the spectrum shown in fig (|5|). 




Figure 7: The cooling of the universe for the case where PBH masses are 
distributed according to the spectrum displayed in fig (|5|). The epoch of 
EWPT is marked by an arrow, and it takes place at 1.9 x 10~ n sec. 



